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Abstract 


We employ molecular dynamics (MD) simulation and experiment to investigate the 
structure, thermodynamics, and transport of A-methyl-A-butylpyrrolidinium bis(tri- 
fluoromethylsufonyl)imide ([pyrl4][TFSI]), TV-methyl-iV-propylpyrrohdinium bis(fluoro- 
sufonyl)imide ([pyrl3][FSI]), and l-ethyl-3-methylimidazolium boron tetrafluoride ([EM- 
IM][BF 4 ]), as a function of Li-salt mole fraction (0.05 < x Li + < 0.33) and temperature 
(298 K < T < 393 K). Structurally, Li + is shown to be solvated by three anion neigh- 
bors in [pyrl4][TFSI] and four anion neighbors in both [pyrl3][FSI] and [EMIM][BF 4 ], 
and at all levels of x Li + we find the presence of lithium aggregates. Pulsed field gradient 
spin-echo NMR measurements of diffusion and electrochemical impedance spectroscopy 
measurements of ionic conductivity are made for the neat ionic liquids as well as 0.5 
molal solutions of Li-salt in the ionic liquids. Bulk ionic liquid properties (density, 
diffusion, viscosity, and ionic conductivity) are obtained with MD and show excellent 
agreement with experiment. While the diffusion exhibits a systematic decrease with 
increasing x Li +, the contribution of Li + to ionic conductivity increases until reach- 
ing a saturation doping level of x Li + = 0.10. Comparatively, the Li + conductivity of 
[pyrl4][TFSI] is an order of magnitude lower than that of the other liquids, which range 
between 0. 1-0.3 mS/cm. Our transport results also demonstrate the necessity of long 
MD simulation runs (^ 200 ns) required to converge transport properties at room T. 
The differences in Li + transport are reflected in the residence times of Li + with the 
anions (r Ll / _ ), which are revealed to be much larger for [pyrl4] [TFSI] (up to 100 ns 
at the highest doping levels) than in either [EMIM][BF 4 ] or [pyrl3][FSI]. Finally, to 
comment on the relative kinetics of Li + transport in each liquid, we find that while the 
net motion of Li + with its solvation shell (vehicular) significantly contributes to net 
diffusion in all liquids, the importance of transport through anion exchange (hopping) 
increases at high x Li + and in liquids with large anions. 
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Introduction 


Over the past decade, there has been extensive research concerning the development of 
rechargeable secondary batteries with high energy density . 1,2 Currently, lithium ion batter- 
ies with organic electrolytes represent the state of the art, though electrolyte stability and 
a maximum energy density of around 200 W-h/kg remain limiting factors . 3 A potential ad- 
vancement beyond current lithium ion technology is the use of lithium metal anodes, which 
would ideally lead to a 5-fold gain in energy density. However, against these anodes, con- 
ventional organic electrolytes are unstable and lead to both poor cyclability and potentially 
hazardous lithium dendrite formation . 4 As an alternative, room temperature ionic liquids, 
which have low volatility, large electrochemical windows, and the ability to solvate Li-salts, 
have been shown to form stable solid-electrolyte interfaces (SEI) against lithium metal, con- 
currently allowing improved cycling and suppression of problematic dendrites . 5-13 In this 
regard, the most promising ionic liquids generally contain either the bis(trifluoromethane- 
sufonyl)imide ([TFSI]) anion or the bis(fluorosufonyl)imide ([FSI]) anion, both of which have 
been shown to allow stable cycling with lithium metal anodes . 5-12 Comparatively, the sta- 
bility of [TFSI] is superior to [FSI ], 11,14 while the discharge rate of [FSI], ~10 mA/cm 2 , 10,11 
has been reported to be larger than [TFSI], 1-1.75 mA/cm 2 . 9,12,15 High discharge rates and 
reasonable stability can also be obtained with mixtures of ionic liquids having small anions, 
notably boron tetrafluoride ([BF 4]), 12,13,16 with organic additives such as vinylene carbonate 
(VC), though cycling in such cases may be generally limited . 12,16 

Bulk transport properties, including viscosity, diffusion, and ionic conductivity, of the 
ionic liquids are important for electrochemical performance, with low viscosity and high 
diffusion being key to high discharge rates. However, while the transport properties of 
neat ionic liquids have been characterized in detail , 17-22 mixtures of ionic liquids with Li- 
salts have not been as thoroughly investigated, with most studies measuring properties at 
select values of x Li +. 5,6 > 9, 13,14,21 Indeed, due to the many possible x Li + and T conditions, 
providing a complete characterization of the influence of Li-salt is difficult, and only a few 
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liquids have been studied comprehensively in this manner . 23-25 Notably, extensive pulsed- 
field- gradient spin-echo NMR studies reported in the works of Hayamizu et al. 24 and Nicotera 
et al. 25 for l-ethyl-3-methylimidazolium ([EMIM]) [BF 4 ] doped with LifBF^ and N-methyl- 
N-propylpyrrolidinium ([pyrl3]) [TFSI] doped with Li[TFSI], respectively, have mapped the 
diffusion and transference numbers as a function of x Li + . 

Complementary to bulk property measurements, the solvation environment of Li + can 
provide important insights into a given liquid. As an example, high energy X-ray diffraction 
(HEXRD) can be used to obtain the liquid structure factor, which in combination with 
molecular dynamics (MD) simulation is an effective tool for finding intermolecular distances. 
HEXRD measurements have been performed to find the relative intermolecular separations 
for ionic liquids with the [TFSI] and [FSI] anions 26-28 and have, interestingly, been adapted 
to characterize Li clusters . 28 Finer structural aspects, such as the binding configurations 
in the Li + solvation shell, are more difficult to obtain and generally require IR or Raman 
measurements in combination with density functional theory (DFT) computations . 27,29-32 
Such procedures have been successfully applied, again, to [TFSI] and [FSI] to identify the 
conformational equilibrium between their cis and trans rotamers 29-31 as well as to determine 
their solvation numbers and preferential bonding states with Li + , 27,32 which include one- fold 
(monodentate) and two-fold binding (bidentate). 

As a supplement to experimental characterization, MD is an effective tool to evaluate 
many fundamental structural, thermodynamic, and transport properties while allowing an 
atom-level view of the structure. Various force fields have been recently parameterized for the 
explicit purpose of ionic liquid simulation, and many of these have matured to the point of 
producing experimentally quantitative thermodynamic properties . 33-37 In particular, polar- 
izable force fields, where each atom is given an environment-dependent atomic dipole, stand 
out as systematically providing superior transport properties in neat ionic liquids . 34,38,39 
More recently, ionic liquid force fields have been adapted to treat the important case of 
mixtures with Li-salts . 40-43 These force fields have been primarily applied to a range of sys- 
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terns, including Li[TFSI] in liquids with pyrrolidinium or imidazolium cations and [TFSI] 
anions, 41,42,44,45 Li[FSI] in [pyr 14] [FSI] , 45 Li[BF 4 ] in liquids with i mi dazolium cations and 
[BF 4 ] anions, 43 and Li[PF 6 ] in [EMIM][PF 6 ]. 46 Such studies generally focus on structure at 
low x Li +, where clustering is not influential, and on properties at high-T, where the typi- 
cally glassy ionic liquid dynamics become amenable to MD simulation. Polarizable models, 
as in the case of neat ionic liquids, are expected to provide more accurate results for both 
thermodynamics and transport properties for systems with Li salt additions. 41,44,45 

In the present study, we employ both theoretical and experimental methods to inves- 
tigate the influence of Li-salt doping on the properties of three potential ionic liquid elec- 
trolytes. The three ionic liquids of interest (see Fig. 1) are IV-methyl-IV-butylpyrrolidinium 
bis (trifluoromethylsufonyl) imide, iV-methyl- .A-propylpyrrolidinium bis (fluorosufonyl) imide , 
and l-ethyl-3-methylimidazolium boron tetrafluoride, denoted as [pyr 14] [TFSI] , [pyr 13] [FSI], 
and [EMIM][BF 4 ], respectively. These ionic liquids have been selected primarily because of 
their stability against Li-metal anodes. 8,9,12 However, their high ionic conductivity and low 
viscosity 21,22,24 make them promising electrolytes for conventional lithium ion batteries 5,14,47 
as well. 

Theoretically, we perform a comprehensive set of polarizable MD simulations, using the 
Atomistic Polarizable Potential for Liquids, Electrolytes, & Polymers (APPLE&P), 33,34,40 to 
investigate the influence of temperature and lithium doping on the structure, thermodynam- 
ics, and transport properties of these liquids. Li-salt is introduced into each ionic liquid using 
the host anion (i.e [pyrl4][TFSI] would have the Li[TFSI] salt), and we probe the influence 
of Li-salt doping from small molar fractions (0.05) to the glassy upper limit of liquid stability 
(~0.33). Importantly, we analyze the ionic liquid properties from room temperature up to 
120 °C at all levels of Li-doping, and, in this regard, our study requires very long simulation 
times (~200 ns) to accurately determine the transport coefficients in low-T ionic liquids. 

In addition, we have performed a number of room temperature experiments on both the 
neat ionic liquids and 0.5 molal solutions of Li-salt in the ionic liquids. Our experiments 
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include measures of density, pulsed field gradient spin-echo NMR measurements of diffusion 
coefficients, and electrochemical impedance spectroscopy measurements of ionic conductivity. 


Methods 

APPLE&P Interatomic Potential 

To accurately describe the complex, environment-dependent electrostatic interactions of 
ionic liquids with MD simulation, we employ the APPLE&P force field as parameterized 
by Borodin and coworkers. 33:3110 The energy in APPLE&P is decomposed into intramolecu- 
lar, repulsive-dispersive, and electrostatic terms. The intramolecular interactions consist of 
harmonic bonds, angles, and improper dihedrals, while torsional dihedrals are represented by 
an order - N harmonic series. The repulsive-dispersive interactions are treated with a mod- 
ified form of the exponential-6 Buckingham potential, which includes a correction term to 
account for close atomic approaches, given by 

u ri ( r tj ) = Aqe-*** - C tj rJ + D tj . (1) 

Furthermore, each atom is assigned a static parital point charge, all of which mutually 
interact through Coulomb’s law. While many of these interaction terms are standard to MD 
simulation, the unique aspect of APPLE&P that allows for a high-fidelity representation of 
ionic liquid energetics is the inclusion of self-consistent atomic polarization. This is achieved 
by associating an atomic polarizability (a) to each atom, which defines the atomic dipole 
moment as proportional to the local electric field, n- = oEfr,). Self-consistency is then 
achieved by updating the electric field with the evolved dipole contributions and recalculating 
the atomic dipole moment until changes in polarization energy, given by 

u^(n) = -\^ • E°(rO, (2) 
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where E°(rj) is the electrostatic field contribution from permanent charges only, reach a 
requisite tolerance. Specific to the [FSI] and [TFSI] anions investigated in this work and in 
addition to the atom center terms, a single lone-pair force center is associated with the S-N-S 
bridge, and the position of this center is instantaneously determined during the simulations 
as a given distance, t = 1.4 A, along a bisector of the S-N-S angle. 

Molecular Dynamics Simulations 

The majority of the MD simulations in the present work are performed with the Large-scale 
Atomic / Molecular Massively Parallel Simulator (LAMMPS). 48 We have modified LAMMPS 
to include the following functionalities required by APPLE&P: lone pairs force centers for 
the [TFSI] and [FSI] anions; self-consistent update of atomic dipoles using the previously 
discussed atomic polarizability; and energy, force, and virial terms that arise from long 
range charge-dipole interactions. For purposes of validation and statistical error estimation, 
a small subset of simulations on [pyrl4][TFSI] are duplicated with Lucretius 49 an simulation 
package specifically tailored to perform MD simulation using the APPLE&P. The simulations 
are integrated using a 3-2-1 reversible reference system propagation algorithm (rRESPA), 50 
which allows the use of a large 3 fs timestep for the long-range Ewald summation of the 
charge and dipole interactions; a 1.5 fs timestep for local charge, repulsive-dispersive, and 
dihedral interactions; and 0.5 fs for the remaining bonds, angles, and improper dihedrals. The 
cutoff used for non-bonded interactions is 12 A for [pyrl4] [TFSI] and [pyrl3] [FSI] systems 
and 11 A for [EMIM][BF 4 ] systems; the cutoff used for non-bonded interactions in the 1.5 
fs rRESPA step is 7 A. For further efficiency, a coupled reaction field/Ewald summation 
formalism is employed in the outer rRESPA partition to evaluate the dipole-dipole (reaction 
field) and dipole-charge (Ewald summation) long range interactions; these computations are 
performed iteratively in the outer partition until the change in the dipole energy reaches a 
convergence criterion of 10 -9 kcal/mol. 

The systems used in our MD simulations contain 144 cation/anion pairs for the [pyrl4]- 
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[TFSI] ionic liquid and 216 pairs for both the [pyrl3][FSI] and [EMIM][BF 4 ] ionic liquids, and 
Li-salt is introduced into each system by replacing a given number of cations with lithium 
ions. At T = 298 K, we perform MD simulations across a range of Li-salt mole fractions 
(x Li + from 0.05 to 0.33), while at higher- T we perform calculations at select values of x Li + 
that correspond to specific experimental conditions. To equilibrate each system, we perform 
simulations at ambient P and the desired T using the Nose-Hoover thermostat and barostat. 
The thermalization is carried out over 20 ns at each T, which we have found to be sufficient 
for density convergence in the system with the lowest temperature and the highest Li-salt 
doping. 

There are additional challenges to the MD simulation of ionic liquids with Li-salt doping, 
primarily due to the glassy nature of the systems at T < 333 K. This is especially relevant 
when considering the calculation of transport coefficients, which are well known to converge 
more slowly than other properties like density. To ensure convergence, we progressively 
compute the transport properties by incrementally adding 30-50 ns of simulation statistics 
to the transport average until the change in the property of interest is less than 5%. In this 
way, we have determined that simulations of 200 ns are required to converge the transport 
properties for the T = 298 K and 333 K systems, while 90-120 ns simulation times appear 
more than sufficient to converge properties at higher T. 

Experimental Techniques 

We have performed extensive experimental characterization of [pyrl4][TFSI], [EMIMjfBF^, 
and [pyrl3][FSI] in both neat samples and those having 0.5 moles of Li-salt per kilogram of 
ionic liquid, denoted as a 0.5 molal solution (m). For these liquids we have measured the 
density, diffusion and ionic conductivity. 

Density is measured under an argon atmosphere (humidity < 1 ppm) in a glovebox at 
room temperature, using a 5 cm 3 graduated cylinder (0.1 cm 3 graduations). Mass is measured 
to the nearest milligram with a calculated uncertainty of <3%. 


The diffusion coefficients of cation and anion species are measured by pulse-field gradient 
spin echo NMR. The signals of the methyl 1 H and 19 F are used to determine the self-diffusion 
coefficients (D + and D“) of the cations and anions, respectively, and in the presence of 
lithium salt, the signal of 7 Li is used to determine the Li + diffusion coefficient. In the NMR 
experiments, D is obtained using 

ln(A/A 0 ) = -D 7 2 (V - S/3)S 2 g 2 , (3) 

where A and A 0 represent the signal integrals in the presence of and absence of the pulsed- 
field gradient, 7 is the nuclear magnetogyric ratio, V is the time interval between the two 
gradient pulses, 6 is the gradient pulse width, and g is the gradient magnitude. 

Samples are loaded into a 5 mm cylindrical thin wall 9” NMR tube inside a glovebox 
and at room temperature, and then the NMR tube was sealed tightly with the cup coming 
with the tube as the tube was purchased. A total of four samples were prepared and loaded 
into the tubes for NMR measurements: [pyrl4][TFSI]; 0.5 m Li[TFSI] in [pyrl4][TFSI]; 
[EMIMjfBF/t]; 0.5 m Li[BF 4 ] in [EMIM][BF 4 ]. All samples were made with pure materials 
and no d-solvent was added. Although the samples inside the NMR tubes were cupped and 
tightly sealed, the samples were also kept in a desiccator to prevent moisture pick-up before 
measurement. 

The samples are inserted into 5 mm BB (broad band) PFG (pulsed field gradient) probe 
equipped with a triple axis gradient amplifier, and all the measurements were collected on 
a Varian VNMRS 500 MHz spectrometer. The grade type “ppp" was used to achieve the 
highest gradient strength of 130 G/cm, which is the result of using two PFG amplifiers 
together with the grade type selection. A doped D 2 0 standard sample from Agilent was 
used to calibrate the gradients before the experiments are carried out. In addition, before 
each experiment, the 90° pulse width was determined on each nuclei ( 1 H, 19 F, and 7 Li) for 
every sample. Typically, 32-64 transients were used to collect the data for each experiment 
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depending on the signal to noise ratio. The diffusion coefficients of Li + , cations, and anions 
are measured at 30° C, with sample temperature being controlled by a variable temperature 
control unit using a heated dry air. The experimental errors in the diffusion coefficient 
measurements are estimated to be less than 3%. 

Ionic conductivities are determined by electrochemical impedance spectroscopy (EIS) 
methods, using a Solarton 1250 frequency response analyzer coupled to a Solartron 1286 
electrochemical interface. The operation of these instruments is automated using Scribner 
and Associates’ Zpolt® software for Windows. Impedance is measured using a 10 mV ampli- 
tude signal (55 kHz to 0.1 Hz) with the cell held at the open-circuit potential. Such methods 
have been previously used to measure the conductivity of polymer electrolytes. 51,52 The con- 
ductivity fixture, displayed in Figure S6 of the supplementary, uses 1-inch diameter, AISI 
304 stainless steel electrodes, mounted in rigid, electrically insulating plates. The electrodes 
are lapped flat with the inside face of the insulating plate, and polished with 600 grit emery 
cloth. Steel dowel pins control the alignment of the electrodes, and inter-electrode spacing 
is established by a slotted, polypropylene mask, which confines the liquid sample between 
the electrodes. The electrode spacing is typically 1 mm. 

Material handling and cell assembly are performed inside an inert atmosphere glovebox 
to minimize the absorption of water (humidity < 1 ppm). After assembly, the fixtures are 
sealed in modified Nalgene™ jars before being transferred out of the glovebox atmosphere 
for testing. Measurements are typically completed between 20°C and 80°C, with temperature 
being allowed to stabilize at ±0.1°C before testing. Electrode spacing, L, and area, A, are 
controlled by the fixture, and electrochemical impedance data provides a measure of sample 
bulk-resistence, Rt,. Ionic conductivity, A, is calculated using A = L/ARi,. The error in the 
calculated conductivity is ±10%, based on tolerances for L, A, and R/,. 
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Structural Implications of Li-doping 


Ionic Liquid Coordination and Conformations 

To assess the influence of Li-salt on the molecular coordination of the ionic liquid, we com- 
pute the center of mass radial distribution function, g(r), between the cations and anions, 
Fig. 2. Interestingly, increased Li-salt doping induces only subtle changes in the cation-cation 
distributions, g ++ (r), and cation-anion distributions, g +- (r), with the relative magnitude of 
the first solvation shell slightly shifting downward. Alternatively, the first solvation shell of 
the anion-anion distributions, g (r), exhibits a significant split into two individual peaks, 
which occurs such that there is a well defined peak at small r and a more diffuse peak 
at larger r. The former of these peaks represents the closely neighbored anions that are 
strongly bound to Li + , with average neighbor distances of 6.5 A, 6 A, and 5 A for [TFSI], 
[FSI], and [BF 4 ], respectively. The latter peak, which is situated at 12 A, 10 A, and 9 A for 
[TFSI], [FSI], and [BF 4 ], respectively, represents the remaining anions coordinated to cations 
and suggests that the average anion-anion neighbor distance around each ionic liquid cation 
increases with Li-salt doping. 

To further elaborate on the coordination of the ions, we compute the average number of 
anions neighboring each cation, given by n catlon (anion) = Anp f g + ~(r)r 2 dr, where p is the 
number density of anions and the integral is taken over the first solvation shell. Of the neat 
liquids, cations in [pyrl4][TFSI] have the fewest coordinated anions, nt pyrl4 l([TFSI]) = 5.6, 
followed by [EMIM][BF 4 ] with n [ EMIM l([BF 4 ]) = 6.7 and [pyrl3] [FSI] with rJ pyrl3 ]([FSI]) = 7. 
While g +_ (r) apparently changes little with increased Li-salt, n catlon (anion), Fig. 3a, shows 
a pronounced increase, with each cation having on average more than one additional coor- 
dinating anion at the highest lithium mole fraction, x Li + = 0.33. This qualitatively agrees 
with the noted splitting behavior of the first solvation shell of g (r), as more anion neigh- 
bors would have increased mutual electrostatic repulsion that would shift their coordinating 
distance to larger r. Of equal importance, this also implies that Li + , on average, maintains 
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fewer neighboring anions than that of an ionic liquid cation, and this duality between cations 
with a large, unfavorable number of anions neighbors and Li + with fewer, strongly bound 
anions could be an energetic penalty to Li-solubility at high levels of doping. 

An additional point of interest is the fact that the larger anions, [TFSI] and [FSI], as- 
sume distinct cis and trans rotamers that can be markers for changes in binding preference 
brought about by Li-salt doping. In the ionic liquids, we have mapped out the distributions 
of the C-S- • -S-C and F-S- • -S-F dihedrals, </, for [TFSI] and [FSI], respectively, noting that 
cis conformations occur when <p approaches zero and trans conformations occur at (f) close 
to 180°; the native distributions at room temperature are mapped out in Fig. Si of the sup- 
plemental material . 53 It appears that the conformational distribution of cis and trans anions 
are slightly influenced by Li-salt levels, shown in Fig. 3b, which presents the conformational 
distribution at T = 298 K as a function of x Li +. While cis- [TFSI] is slightly more probable 
(54 %) and cis-[FSI] is slightly less probable (36 %) than trans conformations at x Li + = 0 , the 
population of cis rotamers increases with increasing Li-salt doping. This agrees with various 
coupled Raman/DFT investigations of these anions, showing that cis and trans conformers 
of [FSI ] 26,54 and [TFSI ] 29,31 exist in equilibrium in liquid phases, with the population of cis 
increasing with Li-salt doping . 30 

Li + Solvation Shell 

As the most general measure of the lithium solvation structure, we compute g(r) between 
Li + and the center of mass of the ionic liquid anions as well as Li + and the atomic species 
that mediate binding with the anion, which is O for [(T)FSI] and F for [BF 4 ]. As shown in 
Fig. 4a-c, [TFSI] and [FSI] show split initial peaks centered around 4 A, while [BF 4 ] has a 
single narrow peak around 3 A that changes little with increasing lithiation. The split first 
solvation shell of [TFSI] and [FSI] is primarily an isomeric effect, which leads to cis-[(T)FSI] 
having a closer center of mass distance to bound Li + than trans rotamers; however, the 
variation in the relative heights of the peaks with x Li + also signify changes in Li + /anion 
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bonding configuration that will be discussed more fully later. The corresponding atomic 
g(r) reveal lithium/atom binding distances close to 2 A, Fig. 4d-e, and further show no 
net change in structure with increasing x Li + , though the heights of the solvation shells shift 
slightly downward. The slight changes in g(r) are again more prominent when considering 
the total number of neighbors, Table 1. As previously speculated, Li + maintains far fewer 
neighbors than ionic liquid cations, with n Ll+ (anion) being less than four in [pyrl4] [TFSI] 
and four or greater for [pyrl3] [FSI] and [EMIMjfBFzt]. Increasing the lithium-doping levels 
generally increases the number of relative anion neighbors, with the strongest effect being 
seen in [pyrl4] [TFSI] and [pyr 13] [FSI] , while the atomic coordination number is generally 
maintained, or barely increasing as with [FSI]. This of course implies that the number of Li + 
bonds per anion decreases at higher x Li +. 

Expanding on this, we have identified the preferential binding states of Li + with the ionic 
liquid anions, representative structures and binding distances of which are given in Fig. S2. 53 
For each of our systems, we find states in the bulk liquid that correspond to Li + having a 
single bond with each anion, or a monodentate configuration ( k 1 ) , as well as those having 
a two- fold bond with each anion, or a bidentate configuration (fid). The binding distances 
increase in magnitude with Li + — Cd FSI ^ > Li + _ 0 [ TFSI 1 > Li + — F^ 4 , and, consequently, 
this is reflected in the binding energetics, which follow Li + -[FSI] < Li + -[TFSI] < Li + -[BF 4 ], 
A more detailed look at the binding energetics reveals fid binding is generally stronger than 
fid states, which are not stable in the gas phase for [TFSI] and [FSI], though solvation effects 
in MD simulations apparently stabilize fid in [TFSI] and [FSI] with respect to gas phase. The 
extent of the stabilization of fid can be seen through comparison of the relative proportion of 
fid binding to fid binding, Fig. 5, which reveals that only [TFSI] has appreciable quantities 
of fid binding (38%) at low x Li +, while |BF 4 | and [FSI] have 95% and 96% fid binding, 
respectively. Increasing lithium concentration further increases the proportion of fid to fid 
in all cases, though this effect is strongest in [TFSI] with a 30% increase in fid binding at 
x Li + = 0.33. 
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As a final note on the lithium solvation shell, we depict in Fig. 6 the most likely coordi- 
nation complexes of Li + with neighboring anions as found at T = 298 K and x Li + = 0.05, 
where the influence of Li + aggregates is minimal as shown in the Li + - Li + radial distribution 
functions in Fig. 7. In [pyrl4][TFSI], roughly 66% of the Li + is coordinated by two k 1 anions 
and one n 2 anion, though complexes having two k 2 anions with one n 1 anion (18%) and 
four anions that are exclusively k 1 (12 %) are also competitive. In the other liquids, Li + is 
almost predominantly coordinated by four k 1 anions, with 80% and 83% of Li + assuming 
this coordination in [pyrl3][FSI] and [EMIM][BF 4 ], respectively. The remaining coordination 
complexes in these two liquids are generally composed of one k 2 anion and three k 1 anions, 
which leads to Li + being atomically 5-coordinated. The T-dependence of the coordination 
complexes is weak, though monodentate configurations become marginally more favorable at 
high-T ; a full analysis of the likelihood of a given coordination complex found in our systems 
as a function of T is further illustrated in Fig. S3. 53 

Qualitatively, the solvation structures presented here agree with available MD simulations 
of similar systems that use both polarizable and non-polarizable force fields. 42,44,45,55 When 
compared to experiment at higher x Li +, however, our predicted lithium solvation structures 
for [pyrl4] [TFSI] and [pyr 13] [FSI] have higher anion coordination numbers. Infrared and 
Raman spectroscopy coupled with DFT for Li[TFSI] and Li[FSI] in imidazolium-based ionic 
liquids 27,32 have been interpreted as showing that [TFSI] prefers two coordination while [FSI] 
prefers three coordination. The experimentally reported coordination numbers are evaluated 
in systems having x Li + > 0.05, which, as shown in the following section, form significant 
quantities of Li + aggregates. It is well known that such aggregates have the net effect of 
reducing the Li + /anion coordination number, 32 and, as it is not simple to spectroscopically 
characterize their contribution to the total measure of coordination, could be the primary 
reason for the discrepancy between our coordination number and that of experiment. 


14 


Li + Clusters 


As previously mentioned, spectroscopic experiments 32 have noted that the Li + /[TFSI] coor- 
dination number decreases at high x Li + , which has been attributed to the tendency of Li + 
to aggregate into clusters - networks of lithium ions bound by shared, or bridging, anions. 
The signature of such clusters has additionally been noted in HEXRD experiments 28 as well 
as MD simulations. 41 Here, too, we observe the presence of Li + clusters and attribute the 
majority of the aforementioned x Li +-mediated structure dependence to their formation. To 
put the magnitude of this effect into perspective, we first look at the probability of finding 
a cluster with a given number of bridged Li + ions, N Li +, as a function of x Li +, shown in 
Fig. 8. In all cases, as x Li + increases, the average cluster size increases, with the probability 
of finding a Li + not associated with a cluster decreasing from ~80-90% at x Li + = 0.10 to 
~ 50% at x Li + = 0.33. There also appears to be a correlation between anion size and cluster 
size, with larger anion sizes leading to progressively fewer, and smaller, clusters. At high 
doping levels the largest clusters observed in our systems are N Li + = 5 in [pyrl4] [TFSI] and 
N Li + = 6 in [pyrl3] [FSI] and [EMIM] [BF 4 ] . 

It is further instructive to consider the radial distribution function of Li + with Li + , which 
are displayed in Fig. 7 for multiple values of x Li +. Of primary importance are the well- 
defined peaks in the initial r < 10 A region that signify Li + aggregation. From inspection of 
simulation trajectories, peaks between ~5-10 A represent lithium ions in small clusters, where 
possible anion binding sites are not saturated and allow configurational entropy to maximize 
Li + separation. As larger clusters form, internal anions become more highly networked until 
all possible binding sites, including those imposing small Li + separations, are saturated, 
which is represented by the peaks in the 4-5 A regions. The initial peak is not present in 
[pyrl4] [TFSI] with x Li + = 0.05, suggesting no large clusters have formed, while the peak is 
present in both [pyrl3] [FSI] and [EMIM][BF 4 ]. In all cases, the relative height of the first 
peak to the second peak increases with x Li + and indicates the formation of large clusters. 

As previously illustrated in Figs. 3a and 5, [TFSI] and [FSI] tend toward k 1 binding and 
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cis-rotamers as lithium-salt doping increases. To examine the contribution of Li + aggrega- 
tion to these changes, we have evaluated the probability of finding the cis rotamer and k 1 
binding in anions coordinated with multiple Li + , the N Li + = 5 cluster represented in Fig. 9 
showing the possible anion coordinations (one Li + being the minimum and four Li + being 
the maximum). For both the [TFSI] and [FSI] anions, as coordination increases, we find 
that the anions tend to primarily assume k 1 binding states and cis configurations. At the 
highest coordination, both anions are 100% k 1 , each available O binding to a Li + , while 86% 
of [TFSI] anions and 82% or [FSI] anions assume the cis conformation. These changes mirror 
experimental studies of the condensed phases of Li[FSI] and Li[TFSI], which also generally 
trend toward monodentate binding 56 and cis-rotamers. 57,58 Further evidence of the influence 
of clusters is provided in Figs. S4 and S5 of the supplementary, which explicitly evaluate 
the cluster contributions to changes in cis-rotamers and k 1 binding in the same manner as 
Figs. 3a and 5. 


Thermodynamic and Transport Properties 

Density 

The density (p) of our ionic liquid systems is shown in Fig. 10 for both the neat systems as 
well as those doped with Li-salt. We find that density agrees to within 1% of experimental 
measures for both the neat and Li-salt doped liquids. Furthermore, the densities of the 
neat ionic liquids show the typical linear decrease with rising T. We find that p trends 
upward with the size of the ionic species, with p[py rl4 H TFSI ] > p[pyri3][FSi] > p[EMiM][BF 4 ]^ w j 1 j c ] 1 
is sensible in light of the fact that average ion separations, on the basis of the g(r) results 
in Fig. 2, do not appreciably change from one liquid to another. The introduction of Li- 
salt, Fig. 10b, results in a near-linear increase up to the highest doping levels investigated, 
x Li + = 0.33. Of the three liquids [pyrl3] [FSI] shows the least linearity and, correspondingly, 
the largest p response to Li-salt doping, with p[py rl3 H FSI l approaching that of [pyrl4] [TFSI] 
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at high x Li +. The relatively large p increase is indicative of the efficient packing of [FSI] 
around [pyrl3] and Li + , which is supported on the basis of Fig. 3a and Table 1, where the 
[pyrl3] [FSI] have the largest values of n catlon (anion) and n Ll+ (anion) at high x Li +. 


Transport Properties 

To compute the shear viscosity (rj xy ) of the ionic liquids, we have employed the Green-Kubo 
relation, 

v r°° 

Vxy = k^fJ 0 dt i ( 4 ) 

where V is system volume, Ub is Boltzmann’s constant, and P xy is shear stress. The viscosity 
of [pyr 14] [TFSI] , [EMIM][BF 4 ], and [pyrl3][FSI] is shown in Fig. 11 for both the neat liquids 
and those having x Li + = 0.10 over the T range 298 K to 393 K. All of our computational 
estimates show agreement to within 20% of available experimental measures. 19,24,45 Addi- 
tionally, rj xy exhibits both a rapid decrease as T increases and a significant increase with 
x Li +. At low-T, [pyr 14] [TFSI] has the largest values of rj xy , which is slightly less than a 
factor of two greater than that [EMIM][BF 4 ] and [pyrl3][FSI], which has the lowest tj xy . At 
high-T, r\ xy of all systems are comparable, though, interestingly, at this point Li-salt doping 
has less influence on q xy of [pyrl4] [TFSI] and [pyr 13] [FSI] than that of [EMIM][BF 4 ]. 

We compute diffusion coefficients with the Einstein relation, 


= lim 


1 d 


t-Hx> 6n^ dt 


i r <(*) 


, £=1 



(5) 


where the diffusion coefficient, D, for molecular species fd is estimated using the position, r, 
of each of the iF atoms belonging to species fd. In all of our quoted values of D, we have 
applied a hydrodynamic correction, which accounts for the influence of the finite simulation 
cell on the diffusion coefficient. The correction, D FSC , is calculated according to 


j^fsc _ 2.837fceT 
67 t r/L ’ 


(6) 
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with L being the average MD simulation cell length. Measures for the diffusion coefficients 
of [pyrl4][TFSI], [EMIM][BF 4 ], and [pyrl3][FSI], both neat and with x Li + = 0.10, are given 
in Fig. 12 for the the T range 298 K to 393 K. Our computational values of D for the 
neat ionic liquid cations and anions in [EMIM][BF 4 ] match available PFG-NMR measures to 
within 5-15%, 24,45 while variations in experimental measures of D in neat [pyrl4][TFSI] lead 
to broad range of potential error (5-50 %). 19,21,45 In general, our estimates of Li + diffusion 
are also within this error tolerance as well; however, we note that at high-T, D of Li + in 
the [pyrl4][TFSI] system seems to be underestimated by a factor of ~2-3, with different 
experiments showing large variations and predicting the Li + diffusion to rapidly approach 
that of [TFSI] at T > 333 K. For [pyrl4] [TFSI] , Fig. 12a and b, Li-salt doping reduces the 
diffusion coefficients of all molecular species, with the general relative magnitudes of diffusion 
following D Ll+ < D' tfsi I < D^ pyrl4 1 + . The same general trends apply to [EMIM][BF 4 ], 
Fig. 12c and d, as well as [pyrl3][FSI], Fig. 12e and f, and the relative magnitudes of diffusion 
follow [EMIM][BF 4 ] > [pyrl3][FSI] > [pyr 14] [TFSI] . Interestingly, in contrast to the other 
liquids, the cations in [pyrl3][FSI] tend to have a lower diffusion than the anions. 

We compute ionic conduction (A) using an expression akin to the Einstein relation for D, 


A = 


lim — — — — 

t-s-oo QVknT 




(7) 


where e is the elementary unit of charge and the sum is taken over all atoms. For [pyrl4] [TFSI] 
and [EMIM][BF 4 ], Fig. 13 a and b, our computed A generally shows agreement to within 15% 
of experiment, 19,24,53 with the major exception being T = 298 K, where the A of [pyrl4][TFSI] 
is a lower bound with an error ~40%. For [pyr 13] [FSI] , Fig. 13 c, on the other hand, there is 
considerable scatter in the experimental data, our computed values agreeing to within 10% 
of the lowest estimates and to within 30% of the highest estimates. The significant scatter in 
[pyrl3] [FSI] A measurements has been previously noted by Zhou and coworkers, 22 who sug- 
gest moisture contamination as a possible cause of higher values of A. Additionally, for Li-salt 
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doped systems, underestimates of A could be a repercussion of the slight overbidding of the 
force field for [TFSI] and [FSI] to Li + when compared to quantum chemical computations. 59 
As expected, T has a strong influence on A, which increases by an order of magnitude from 
298 K to 393 K, and the addition of Li-salt decreases the net ionic conduction. At high-T, 
the ionic conduction of [EMIM][BF 4 ] is a factor ~4 greater than [pyrl4][TFS] and a factor 
of ~ 2 greater than [pyrl3] [FSI] , while this increases to ~ 10 for [pyrl4] [TFSI] at lower T. 
Additional information on the experimental measurement of A is presented in Figure S7 of 
the supplementary material. 

In Fig. 14, we compare our experimental values for A of the ionic liquids (0.5 m solutions 
of Li-salt in the ionic liquids) against that of a standard Li-ion electrolyte (1 mole of Li[PF 6 ] 
per liter of a 1:1 ethylene carbonate and dimethyl carbonate mixture) to asses the suitability 
of our ionic liquids for batteries. The temperature dependence of ionic conductivity is fitted 
to a form of the Vogel- Tamman-Fulcher relationship, 

= w exp ’ <8) 

where A, E a , and T a are fitted parameters and R is the ideal gas constant. The conductivity 
of 0.5 m Li[BF 4 ] in [EMIM][BF 4 ] and 0.5 m Li[FSI] in [pyrl3] [FSI] compare favorably with the 
Li-ion electrolyte at room temperature, but both fall off more rapidly at reduced temperature. 
A of 0.5 m Li[TFSI] in [pyrl4] [TFSI] is nearly an order of magnitude lower than that of the 
other ionic liquids, as well as that of the standard electrolyte, at low and moderate values 
of T. It has, however, been suggested that additives might be a possible route to tune ionic 
conductivity values of potential ionic liquid electrolytes. 

We additionally compute room temperature properties especially important for electro- 
chemical applications, which include r] xy , the diffusion coefficient of Li + (D Ll+ ), and the Li + 
component of ionic conduction (A Ll+ ), for x Li + in the range 0.05 to 0.33, Fig. 15. While ex- 
pressions for r] xy and D Ll+ have been provided, we estimate A Ll+ from a ratio of the diffusion 
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coefficients and the net ionic conductivity, 


A Li+ 


= A 


N Li+ D Li+ 

N + D + + N“D“ + N Li+ D Li+ 


(9) 


where N + , N _ , and N Ll+ are the number of cations, anions, and Li + ions, respectively, in 
our simulation cell. The errors in our computed values of r ] xy , D Ll+ , and A Ll+ generally range 
between 10-15%, though our computed value of A Ll+ at x Li + for [pyrl4][TFSI] is off by a 
factor of two, which is a propagation of the previously mentioned underestimation of A in 
this system at low-T. We note a few interesting trends with respect to increasing Li-doping, 
including increasing rj xy , decreasing D Ll+ , and a plateau in A Ll+ . For each liquid, rj xy and 
D Ll+ increases and decreases, respectively, by an order of magnitude as x Li + rises from 0.05 
to 0.33. At x Li + = 0.33 the [pyrl4] [TFSI] system has the largest viscosity at f448 cP, ~4 
times larger than the other systems, and lowest D Ll+ at 0.006xl0 _1 °, ~4 times smaller than 
the other systems. While the increase and decrease in ri xy and D Ll+ , respectively, are sensible 
and likely correspond to the increase in density with x Li + noted in Fig. fO, the plateau in 
A Li+ can be attributed to a balance between the diffusion of lithium, which decreases with 
higher x Li +, and the total number N Li +, which increases. In terms of A Ll+ , [pyrI4][TFSI] 
again has the poorest performance at 0.03 mS/cm, which is nearly an order of magnitude 
lower than A Ll+ for [EMIM][BF 4 ] at 0.26 mS/cm and [pyrl3][FSI] at 0.17 mS/cm. 

To better understand the error tolerance of our low-T computational measures of trans- 
port, we now compare our transport values in Figs. 11, 12, and 15 to those from additional 
simulations of [pyrl4][TFSI] having x Li + = 0.1 and 0.25 as performed with Lucretius , 49 The 
Lucretius simulations, which are denoted as “luc,” use initial configurations uncorrelated to 
our other simulations and contain 160 pairs. As shown in these figures, D Ll+ , viscosity and 
A Ll+ from the luc simulations are in good agreement with our standard results derived from 
LAMMPS, with the largest deviations being in the diffusion of Li + and viscosity at low-T 
(typical deviations being ~ 10-15 %). The simulation statistics, then, so obtained by the 
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two simulation packages are congruent. The variability in transport is likely less dependent 
on initial configuration than errors related to the length of simulation time, which is further 
discussed in Fig. S8 of the supplementary material on the basis of deviations in low-T D Ll+ 
over the course of hundreds of nanoseconds of statistical averaging. Thus, we strongly sug- 
gest that care be taken to ensure that ample statistics are taken over long simulation times 
to produce converged transport properties in Li-doped room-T ionic liquids, with hundreds 
of nanoseconds of simulation statistics being necessary to prevent large errors in transport 
estimates. 


Kinetics of Lithium Transport 


It is worthwhile to further analyze the underlying mechanisms that govern the transport, in 
particular that of Li + , in our ionic liquid systems. As a first approach, we compare the A of 
our liquids, as computed from Eq. 7, with an ideal measure of ionic conduction (X uc ) that 
does not account for the correlated motion of the ions, 


X uc 


k B VT 


js ^ions 


(^) 2 D^, 

0=1 


( 10 ) 


where z 13 is the net charge of each ionic species 3. Expressing the effect of correlated ion 
motion as the ratio a = X/X uc , we provide in Table 2 measures of a for the ionic liquids 
at T = 298 K and selected values of x Li +. We find that our estimated a is on the order 
of both experimental measurements, which have yielded 0.39-0.86 for [pyrl3] [TFSI] with 
Li[TFSI] 25,60 and 0.7-0.85 for [EMIM][BF 4 ] with Li[BF 4 ] , 24 as well as alternate MD studies 
of Li-doped ionic liquids. 20,41,44,45 With increased x Li +, we observe a slight increase in a that 
suggests less correlated ion diffusion. 

In order to further probe possible mechanisms of transport, we look to the average resi- 
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dence time of neighboring ions. This is accomplished by defining a neighbor function (H s ' r ) 


H^(R ij,t) 


jo, if Rij > R^ AX| 
(l, if Rij < rJ^ ax| 


( 11 ) 


which determines whether the ions i and j of species type /3 and 7, respectively, are neighbors 
based on their center of mass separation, R ; ; ? , and the average maximum extent of the first 
solvation shell, Rj^ AX ' , which can be obtained from the g(r). The residence time, then, can 
be evaluated from the autocorrelation of H^ 7 , 


„ / Efi Efi h /37 (r ij, t) h /37 (%, 0) \ 

S fr(t) = ( ™ ^ J=1 > , (12) 

\ Ef=i E^i H /37 (%, 0)H /?7 (Rjj, 0) / 

which exhibits a roughly exponential decay, the average decay time of which is r^ 7 between 
ion species j3 and 7. 

The residence time of Li + with anions in each liquid is mapped out in Fig. 16a, which re- 
veals exceptionally long residence times that follow the trend t Li/ 1 tfsi 1 > t Li /I BF/| ] > r Ll /[ FSI ] ^ 
While r Ll / _ in every liquid increases with increased x Li +, the relative rate of increase is 
greatest for [pyrl4][TFSI], with r L, /R FST ] at x Li + = 0.33 reaching 100 ns - fully an order of 
magnitude larger than that of the other liquids. Though a relationship between the two has 
not been previously established, it appears that the increase in r between Li + and the anions 
corresponds to a decrease in correlated motion, as given by increasing a. Additionally, T 
has a marked effect on t L i/ E shown in Fig. 16b for x Li + = 0.10. r L, A decreases by roughly 
an order of magnitude as T is increased from 298 K to 393 K, though we note that r Ll// [ FSI l 
and r Ll /I BF k appear to be converging at 393 K, while r Ll /[ TFSI l i s still a factor of two larger. 

Intuitively, a longer r Ll//_ at higher x Li + would suggest a greater correlated motion; how- 
ever, this would only be the case if ion diffusion rates remained unchanged, which the sup- 
pression of D in Fig. 14 clearly contradicts. Alternatively, the suppressed diffusion and long 
residence time, along with increasing a, could indicate a change in transport mechanism 
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to one in which diffusion occurs with less correlated motion. Consequently, the primary 
mechanisms proposed for Li + diffusion include the highly correlated vehicular mechanism, 
where Li + moves with its solvation shell, and the hopping mechanism, where Li + is pulled 
through the liquid via anion exchange. To compare these mechanisms, we compute the av- 
erage displacement of Li + before an anion exchange event occurs at time r Ll / _ as well as the 
ratio of Li + diffusion before anion exchange (D^) to D Ll+ ; a small ratio of to D Ll+ 
implies more importance of anion exchange on diffusion. In Table 3, we provide D[///D t ' i+ 
at selected conditions along with the Li + displacement measured in numbers of average an- 
ion radii, N^. With increasing x Li +, we see a shift to lower D{^/D Ll+ as well as a lower 
N ' R > , indicating that anion exchange contributes more greatly to net Li + diffusion and that 
Li + is displaced on average less before an anion exchange event, supporting the increase of 
a in Table 2. Based on a comparison of the magnitudes of D^/D Ll+ , the importance of 
the vehicular mechanism follows [EMIM][BF 4 ] > [pyrl3] [FSI] > [pyrl4][TFSI], which also 
implies anion exchange becomes more important as anion size increases. 

As a final evaluation of transport in Li-doped ionic liquids, we look at the influence of 
Li + clusters on both r Ll / _ and D Ll+ . As shown in Table 4, we compute t L i / _ and D Ll+ for 
single Li + ions (denoted clusters of one, N Li += 1) and cluster Li + with common bridging 
anions (N Li + > 1); the average duration of the clusters is also quoted on the basis of the 
Li-Li residence time (r Ll ' /Ll ). We initially note that the magnitude of r Ll/,Ll with respect to 
ionic liquid follows the order previously established for t L i/ ~ and that in all cases t L i / _ i s 
larger for N Li + > 1 systems than N Li + = 1 systems. The magnitude of this effect, however, 
varies markedly by system, with t L i / _ for N Li + > 1 clusters being a factor of ~ 2 and 
~ 3/2 larger than r Ll /“ for N Li + = 1 in [pyrl4][TFSI] and [pyrl3] [FSI] , respectively, while 
both Li + - Anion residence times are comparable for [EMIM][BF 4 ]. Associated with this 
increase in cluster residence time is an increase by a factor of ~ 1.5 in D Ll+ for N Li + = 1 
when compared to N Li + > 1 systems in [pyrl4][TFSI] and [pyrl3] [FSI] , with [EMIM][BF 4 ] 
again being the exception with comparable D Ll+ independent of cluster size. Upon closer 
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inspection of the structures associated with these values, it is noted that the increase in r Ll /“ 
for N Li + > 1 systems can almost completely be attributed to the long-lived bridging anions 
(given by coordination > 1 in Fig. 9), while anions on the exterior of the cluster (given by 
1-coordination in Fig. 9) are exchanged at rates comparable to those associated with N Li + = 
1 clusters. 


Conclusions 

We have presented an comprehensive investigation of the structure, thermodynamics and 
transport properties of the ionic liquids [pyrl4][TFSI], [pyrl3][FSI], and [EMIM][BF 4 ] as 
a function of Li-salt doping and T. Molecular dynamics simulation are performed using 
the polarizable force held, APPLE&P. An analysis of the influence of x Li + on ionic liquid 
structure reveals a complex binding environment in which Li + prefers 4-coordination, which 
is generally mediated through 4 anions having k 1 binding, except for the case of [pyrl4][TFSI] 
which prefers a 3 anion coordination consisting of 2 k, 1 bonds and 2 n 2 bonds. A result of 
adding Li + past a threshold value of x Li + = 0.05 is the mutual coordination of lithium ions to 
the same anion, resulting in lithium clusters. The clustering effect has a significant impact 
on both the structure and transport properties of our ionic liquids, reducing the prevalence 
of n 2 binding in [pyrl4][TFSI] and [pyrl3][FSI] as well as slowing the dynamics of the liquids 
through a substantial increase in density. 

We have evaluated the expected performance of each electrolyte based on their viscosity, 
diffusion, and ionic conduction as well as the individual contribution of Li + to the ionic con- 
duction. MD simulation predictions were found to be in good agreement with pulsed held 
gradient spin-echo NMR measurements of diffusion and electrochemical impedance spec- 
troscopy measurements of ionic conductivity for the neat ionic liquids and 0.5 m solutions 
of Li-salt in the ionic liquids. With increasing Li-salt doping, all liquids exhibit the same 
general trends of increasing viscosity and decreasing diffusion and ionic conductivity. The 
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influence of Li-salt doping is notable, with our T = 298 K systems showing a net increase of 
r] xy and a net decrease of D Ll+ of an order of magnitude as x Li + increases from 0.05 to 0.33. 
Contrasting this behavior, increasing Li-salt doping causes A Ll+ to saturate after a threshold 
of x Li + = 0.15, which is shown to be the combined influence of decreasing diffusion and 
increasing number of Li + . Comparing the properties of our potential electrolytes across the 
entire range of x Li +, a correlation between anion size and transport emerges, with smaller 
anions having smaller viscosity and higher rates of Li + transport. The liquid with the largest 
anion, [pyrl4][TFSI], has the highest viscosity, being a factor of ~4 larger than the other 
liquids, while D Ll+ is a factor of ^4 lower. Without question, [pyr 14] [TFSI] exhibits the 
lowest Li + contribution to ionic conduction, A Ll+ being roughly an order of magnitude lower 
than that of the other liquids, and [EMIMjfBFzt] has the highest A Ll+ , reaching values of 0.30 
mS/cm at high x Li +, though [pyrl3][FSI] is a close lower bound at 0.17 mS/cm. It should 
additionally be noted that our room-T simulation results of Li-doped ionic liquids demon- 
strate the necessity of long MD simulation runs (~ 200 ns) to obtain accurate measures of 
transport. 

An analysis of the kinetics of Li + transport has been performed on the basis of resi- 
dence times, t L i / - , of anions with Li + . At T = 298 K, Li[TFSI] binding persists up to 
100 ns at x Li + =0.33. Li[FSI] and LifBF^ binding, on the other hand, have much smaller 
values of t u '~ at 8 ns and 17 ns, respectively. Increasing T reduces the residence times 
dramatically, exhibiting an order of magnitude decrease from 298 K to 393 K. Using r Ll / _ 
along with diffusion data, we have additionally investigated the competition between the 
hopping mechanism and the vehicular mechanism in the transport of Li + . The contribution 
of Li + /anion exchange to Li + diffusion is estimated to be less that 40 % and increases with 
both increasing x Li + and anion size ([TFSI] > |FSI| > [BF4]). 

Finally, we briefly note that Li + transport is governed by the correlated influence of anion 
size, anion binding strength, and liquid density. For instance, [pyr 14] [TFSI] and [pyrl3][FSI] 
are similar ionic liquids in terms of structure, though the residence time of Li + with [TFSI] 
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is roughly an order of magnitude larger, which inevitably leads to slower diffusion of Li + . 
It is attractive to attribute this to the bonding energy, which is stronger for [TFSI] than 
[FSI] ; however, bonding alone does not completely govern r Ll// , as [BF4] has the strongest 
binding energy with a r Ll /“ several times smaller than that of [TFSI]. What, then, emerges 
is a picture where strong anion binding, large anion size, and high density create a prolonged 
Li[TFSI] solvation structure that leads to slow transport. While LifBFzt] has strong binding, 
the small anion size and low density allow rapid diffusion of Li + with its solvation shell as 
well as a high rate of anion exchange. Finally, the weak binding of Li[FSI] and intermediate 
density of [pyrl3] [FSI] seem to allow a greater rate of anion exchange and result in higher 
Li + mobility. 
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Table 1: Anions and atoms coordinated to Li + at T = 298 K. 


[pyr!4] [TFSI] [pyr!3][FSI] [EMIM][BF 4 ] 


X Li + 

n Li+ ([TFSI]”) 

n Li+ (0) 

n Li+ ([FSI]-) 

n Li+ (0) 

n Li+ ([BF 4 ]-) 

n Li+ (F) 

0.05 

3.3 

4.1 

4.0 

4.1 

4.0 

4.2 

0.15 

3.5 

4.1 

4.1 

4.2 

4.0 

4.2 

0.33 

3.6 

4.1 

4.4 

4.3 

4.1 

4.2 
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Table 2: Ratio of ionic conductivity to uncorrelated ionic conductivity (a) at 

T = 298 K. 


x Li+ 

Q [pyrl4|[TFSI] 

Q/[pyrl3][FSI] 

a [EMIM][BF 4 J 

0.05 

0.64 

0.57 

0.71 

0.15 

0.72 

0.64 

0.78 

0.33 

0.85 

0.73 

0.82 
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Table 3: Ratio of the T = 298 K diffusion of Li + before anion exchange (D^x) to 
the total diffusion of Li + (D Ll+ ). The distance that Li + travels before exchanging 
anions, as measured in numbers of ionic liquid anion radii (N^), is given as 
reference. 


[pyrl4] [TFSI] [pyrl3][FSI] [EMIM] [BF 4 ] 


*Li+ 

d£/d l '* 

n( r > 

Df£/D Li+ 

n< r > 

Df£/D Li+ 

n< r > 

0.05 

0.69 

4.4 

0.81 

3.7 

0.89 

6.1 

0.10 

0.66 

4.2 

0.85 

2.4 

1.07 

5.8 

0.33 

0.59 

3.5 

0.73 

2.0 

0.91 

3.9 
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Table 4: Determination of the influence of clustering on Li + -Anion residence 

time (r Ll//_ ) and Li + diffusion (D Ll+ ) for single Li-ions (N Li + = 1) and those 
belonging to clusters (N Li + > 1) in ionic liquids having x Li + = 0.1 at T = 298 K. 



r Li/Li 

(ns) 

r Li/_ (N Li + = 1) 
(ns) 

r Li /-(N Li+ > 1) 
(ns) 

D Li+ (N Li+ = 1) 
(10 _10 m 2 /s) 

D Li+ (N Li+ > 1) 
(10 _10 m 2 /s) 

[pyrl4] [TFSI] 

81.3 

18.3 

41.1 

0.034 

0.024 

[pyrl3] [FSI] 

8.3 

4.6 

6.8 

0.083 

0.053 

[EMIM][BF 4 ] 

27.7 

11.6 

12.0 

0.078 

0.081 
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(a) 


(b) 




[pyrl 3] + N-methyl-N-proylpyrrolidinium + 

w y 

[EMIM] + 1 -methyl-3-ethylimidazolium + 


trons cis X 

[TFSI] - bis(trifluoromethylsufonyl)imide 

trans cis 

[FSI]" bis(fluorosufonyl)imide 

4 ^ 


[BFJ- boron tetrafluoride 


Figure 1: Representation of the cations (a) and anions (b) of the three ionic liquids consid- 
ered in this work: [pyr!4] [TFSI] , [pyr!3] [FSI] , and [EMIM][BF 4 ]. 
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- [pyrl 4][TFSI] 

— 0.85[pyr14][TFSI] + 0.15Li[TFSI] 
0.67[pyr14][TFSI] + 0.33Li[TFSI] 


- [pyrl 3][FSI] 

— 0.85[pyr13][FSI] +0.15Li[FSI] 
0.67[pyr13][FSI] + 0.33Li[FSI] 


— [EMIM][BF J 

-- 0.85[EMIM][BF 4 ] + 0.15Li[BF 4 ] 

■ - 0.67[EMIM][BFJ + 0.33Li[BF 4 ] 



Figure 2: T — 298 K radial distribution functions for ionic liquid ions as a function of Li-salt 
doping. Distributions are provided for cations with other cations (g ++ ), cations with anions 
(g +_ ), and anions with anions (g ) for (a,d,e) [pyrl4][TFSI], (b,e,h) [pyrl3][FSI], and (c,f,i) 

[EMIM][BF 4 ]. 
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x u+ (mole fraction Li-salt) 


Figure 3: (a) Average number of anion neighbors to each ionic liquid cation, n catlon (anion), 
and (b) probability of cis and trans conformations of anions in the [pyrl4][TFSI] and 
[pyr!3] [FSI] systems a function of Li-salt doping. 


— 0.95[pyr14][TFSI] + 0.05Li[TFSI] 
-- 0.80[pyr14][TFSI] + 0.20Li[TFSI] 

- 0.67[pyr14][TFSI] + 0.33Li[TFSI] 


— 0.95[pyr13][FSI] + 0.05Li[FSI] 
-- 0.80[pyr13][FSI] + 0.20Li[FSI] 

- 0.67[pyr13][FSI] + 0.33Li[FSI] 


— 0.95[EMIM][BFJ + 0.05Li[BFJ 
-- 0.80[EMIM][BFJ + 0.20Li[BFJ 
-- 0.67[EMIM][BFJ + 0.33Li[BF 4 ] 



Figure 4: Radial distributions functions of (a-c) Li + with anions as well as (d-e) Li + with 
atomic binding sites on the respective anions at T = 298 K and as a function of x Li + for our 
ionic liquid systems. 
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Figure 5: Proportion of Li + ligand anions having monodentate bonds as a function of Li-salt 
doping at T = 298 K. 
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Figure 6: Most likely Li + coordination complexes found in (a) [pyrl4][TFSI], (b) 

[pyr!3][FSI], and (c) [EMIM][BF 4 ] at T = 298 K and x Li + = 0.05. 
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5 



r (A) r (A) r (A) 

Figure 7: Radial distribution function of Li + ions with other Li + ions for [pyr 14] [TFSI] , 

[pyrl3][FSI], and [EMIM][BF 4 ] at T = 298 K. The evolution of a large peak at 4 A indicates 
formation of clusters. 
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Figure 8: Li-cluster size statistics given as a function of cluster size at T = 298 K for (a) 

[pyr!4] [TFSI], (b) [pyr!3][FSI], and (c) [EMIM] [BF 4 ], 
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Figure 9: Representative configuration of a lithium cluster taken from a T — 298 K 

[pyr!4][TFSI] system having x Li + =0.33 (F atoms removed for clarity). 
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• [pyr14][TFSI] 

■ [pyrl 3][FSI] 

♦ [EMIM][BFJ 

o [pyrl 4][TFSI] (exp) 
□ [pyrl 3][FSI] (exp) 
o [EMIM][BFJ (exp) 


300 320 340 360 380 400 

T (K) 


• [pyrl 4][TFSI] 

■ [pyrl 3][FSI] 

♦ [EMIM][BFJ 

o [pyr14][TFSI] (exp) 
o [EMIM][BFJ (exp) 


0 0.1 0.2 0.3 0.4 

x Lj * (mol fraction Li-salt) 

Figure 10: Density of (a) neat ionic liquids as a function of system T and (b) ionic liquids 
at T = 298 K as a function of Li-salt doping. MD simulations (filled symbols) are compared 
to available experimental measures(open symbols). 19,26,61 
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1 000/7 (K 1 ) 


Figure 11: Shear viscosity of (a) [pyr 14] [TFSI] , (b) [EMIM][BF 4 ], and (c) [pyrl3][FSI] for 
neat systems and those having x Li + = 0.10 as a function of T. Experimental comparisons 
are included for both systems. 19,22,24,45 
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• [pyrl 4] + 

■ [TFSI] 

A Li + 
a Li + (luc) 
o [pyrl 4] + (exp) 
□ [TFSI]' (exp) 
A Li + (exp) 


• [EMIM] + 

■ [BFJ- 
A Li + 

o [EMIM] + (exp) 
□ [BF 4 ]' (exp) 

A Li + (exp) 


C/3 H 

CM 10 


X 10 
o 



• [pyrl 3] + 
■ [FSI]' 

A Li + 


Figure 12: Diffusion coefficients for [pyrl4][TFSI] in both (a) the neat form and (b) that 

having x Li + = 0.10, [EMIM][BF, 4 ] in both (c) the neat form and (d) that having x Li + = 0.10, 
and [pyrl3] [FSI] in both (e) the neat form and (f) that having x Li + = 0.10 as a function of 
T. MD simulation results (solid symbols) are compared with available experiments (outlined 
symbols). 19,21,24,45 
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Figure 13: Ionic conduction of (a) [pyrl4][TFSI] in both the neat form and that having 

x Li+ = 0.20, (b) [EMIM][BF 4 ] in both the neat form and that having x Li + = 0.10, and (c) 
[pyrl3][FSI] in both the neat form and that having x Li + = 0.15. MD simulation results (solid 
symbols) are compared to experimental measures (outlined symbols). 5 ’ 19,22,24 ’ 53 
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Figure 14: Comparison of experimental measures of A for 0.5 m LifBF^ in [EMIM][BF 4 ], 

0.5 m Li[FSI] in [pyrl3][FSI], and 0.5 m Li[TFSI] in [pyrl4][TFSI]. Experimental measures of 
A for 1 M Li[PF 6 ] in 1:1 ethylene carbonate:dimethyl carbonate (EC:DMC), a typical Li-ion 
electrolyte, are also included. 


49 


© • [pyr14][TFSI] 
o [pyr14][TFSI] (exp) 



0 0.1 0.2 0.3 0.4 

x u+ (mol fraction) 


▲ [EMIM][BFJ 
a [EMIM][BFJ (exp) 


( h ) £_.aaa-- 


0.1 0.2 0.3 

x u+ (mol fraction) 


r (c) 



■ [pyr13][FSI] 

□ [pyr13][FSI] (exp) 


-flr" □ 


0.4 0 


0.1 0.2 0.3 0.4 

x u+ (mol fraction) 


Figure 15: Li-doping dependence of the T = 298 K transport properties, (a-c) shear viscosity, 
(d-f) lithium diffusion, and (g-i) the lithium contribution to ionic conduction, for the three 
ionic liquids considered in this work. MD simulation results (solid symbols) are compared 
to available experiments (outlined symbols). 5 ’ 10,14 ’ 19,21 ’ 24 ’ 53 
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Figure 16: Residence time (Y Ll - / ) of Li + with ionic liquid anions, given (a) as a function of 
x Li + at T = 298 K and (b) as a function of T at x Li + = 0.10. 


51 


SUPPLEMENTARY MATERIAL 


Contents 


1. Structural Analysis 1 

2. Computational and Experimental Thermodynamics 5 

References 9 


1. Structural Analysis 



Figure SI. Distribution of the C-S-S-C dihedral and the F-S-S-F dihe- 
dral in [TFSI] - and [FSI] - , respectively, at T — 298 K. Values of </> close 
to 180° represent trans conformations, while smaller values represent cis 
conformations. 
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(a) 



r Li+ -° = 2.00 A 


(d) 



r Li+ -° = 2.10 A 


Monodentate (k 1 ) 

(b) (c) 

YY* 

r u + -°= 2.03 A 
Bidentate (k 2 ) 



r Li+ -°= 2.21 A 



r Li+ - F = 1.91 A 



r u+ -F= 2 .i2 A 


Figure S2. Depiction of monodentate (a-c) and bidentate (d-f) binding 
states for Li[TFSI], Li[FSI], and Li[BF 4 ], Average bond distances at T = 
298 K are included. 
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Temperature (K) 

Figure S3. Most likely Li + -anion coordination shells found in (a) 
[pyrl4][TFSI], (b) [pyrl3][FSI], and (c) [EMIM][BF 4 ] for a T range 298 
K to 393 K. The values are averaged over simulations at x Li + = 0.05, and 
the shells are defined by the number of monodentate (ft 1 ) and bidentate 
(ft 2 ) anions in the shell. 
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Figure S4. Distribution of cis and trans rotamers as a function of Li-salt 
mole fraction, x Li +, for anions not bonded with Li + (dashed lines) and all 
anions in the system (solid lines). 



Figure S5. Distribution of monodent ate bonds to Li + as a function of Li- 
salt mole fraction, x Li + , for Li ions not associated with a cluster structures 
(dashed lines) and all Li ions in the system (solid lines). 
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2. Computational and Experimental Thermodynamics 

Measurements for neat [EMIM][BF 4 ] exhibit a significant hysteresis in conductivity, 
below -15°C (see Fig. S7). Similar hysteresis has been reported by J. Vila et al. for 
a number of [EMIM]-based ionic liquids. [1] The observed drop in conductivity occurs 
well above the solidification temperature tabulated by Galinski for [EMIM][BF 4 ]. [2] Ionic 
liquids based on [pyr!4][TFSI] were not studied below room temperature in this work. 


stainless-steel 



Figure S6. Conductivity fixture. 


Table SI. Density and diffusion coefficient of [pyrl4] [TFSI] and 
[EMIM][BF 4 ] at T — 298 K for both neat samples and those having 0.5 
m Li-salt doping. 


p (g/cm 3 ) D + (m 2 /s) D (m 2 /s) D Ll (m 2 /s) 


[pyrl4][TFSI] (neat) 

[pyrl4][TFSI] (0.5 m Li-salt) 1.44 
[EMIM][BF 4 ] (neat) 

[EMIM][BF 4 ] (0.5 m Li-salt) 


0.1338E-10 0.1029E-10 
0.0473E-10 0.0270E-10 0.0248E-10 
0.3339E-10 

0.2329E-10 - 0.0765E-10 


1.31 
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1000/T (K' 1 ) 


Figure S7. Conductivity of neat [EMIM][BF 4 ]. 

Table S2. T dependence of the ionic conductivity of [pyrl4] [TFSI] in 
both the neat form and that having 0.5 m Li-salt doping. 

A (mS/cm) 


T (K) 

[pyrl4][TFSI] 

(neat) 

[pyrl4] [TFSI] 
(0.5 m Li-salt) 

293.0 

2.31 

1.22 

313.0 

4.98 

2.85 

333.0 

8.90 

5.55 

353.0 

13.8 

8.54 
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Figure S8. Deviation of computed Li + MSD as a function of initial 
state over the course of a 191 ns simulation at T — 298 K. The difference 
between D Ll taken from the initial trajectory (black solid line) and the 
final trajectory (green dahed-dotted line) is a factor ~ 2. 

Table S3. T dependence of the ionic conductivity of [EMIM][BF 4 ] in 
both the neat form and that having 0.5 m Li-salt doping. 

A (mS/cm) 


T (K) 

[EMIM] [BF 4 ] 
(neat) 

[EMIM][BF 4 ] 
(0.5 m Li-salt) 

352.1 

63.6 

43.6 

334.4 

47.3 

31.1 

332.2 

43.4 

29.7 

314.5 

29.5 

19.3 

313.5 

28.0 

18.4 

297.6 

17.6 

11.0 

295.0 

16.4 

9.7 

294.1 

15.5 

9.6 

275.5 

7.38 

4.0 
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Table S4. T dependence of the ionic conductivity of neat [pyr!3][FSI]. 


T A 


(K) (mS/cm) 

253.7 

1.4 

265.0 

2.6 

274.6 

3.9 

277.3 

4.4 

295.1 

8.0 

295.2 

7.9 

296.5 

8.3 

313.5 

13.1 

318.1 

14.6 

333.1 

20.0 

335.1 

20.1 

353.1 

28.2 


Table S5. T dependence of the ionic conductivity of a 0.5 m solution of 
Li[FSI] in [pyr!3] [FSI]. 


T A 


(K) (mS/cm) 

253.8 

0.9 

255.3 

1.0 

264.8 

1.7 

275.0 

2.9 

284.1 

4.2 

293.8 

6.0 

294.5 

6.2 

294.8 

5.4 

307.3 

8.6 

319.5 

12.7 

333.5 

17.0 

339.3 

19.4 

343.0 

20.7 

353.0 

24.6 
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